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Abstract 

We present numerical results (including full one-loop QCD corrections) for the pro- 
cesses pp and pp — ► W + W~ , W^Z/j* and Z/"f* Z/j* followed by the decay of the 
massive vector bosons into leptons. In addition to their intrinsic importance as tests of 
the standard model, these processes are also backgrounds to conjectured non-standard 
model processes. Because of the small cross sections at the Tevatron, full experimen- 
tal control of these backgrounds will be hard to achieve. This accentuates the need 
for up-to-date theoretical information. A comparison is made with earlier work and 
cross section results are presented for pp collisions at \fs = 2 TeV and pp collisions at 
yfs = 14 TeV. Practical examples of the use of our calculations are presented. 



Submitted to Physical Review D 



1 Introduction 



We present results for the hadronic production of a vector boson pair, including all spin 
correlations in the decay of the final state bosons, qq —>■ V1V2 — ► leptons where Vi = W , Z or 
7*. The calculations are performed in next-to-leading order in as- We implement the helicity 
amplitudes of [Q an d thus extend previous treatments of vector boson pair production (0-0] 
and f|-0]) t° include spin correlations in all the partonic matrix elements. By including the 
decay products in this way it is possible to impose experimental cuts, necessary to compare 
theory with experiment. Some cuts are experimentally necessary and more stringent cuts 
are often useful in order to reduce backgrounds in the search for new physics. Although 
phenomenological predictions including the complete one loop predictions are presented here 
for the first time, this may be a matter of theoretical correctness rather than practical 
importance. The early predictions were performed at 16 and 40 TeV and to a limited extent 
at 1.8 TeV, so an update of the phenomenological results is in any case appropriate. We 
have therefore provided predictions for pp collisions at 2 TeV (Tevatron Run II) and for pp 
collisions at 14 TeV (LHC). Since the early predictions were made, there have been changes 
in a s and the determination of the gluon distribution, especially at small x. We include this 
information by using modern parton distributions. 

Our results are obtained using a Monte Carlo program MCFM which allows the calculation 
of any infra-red finite quantity through order a s . The Monte Carlo program is constructed 
using the method of Ref. || based on the subtraction technique of Ref. ||. We hope to 
provide further details in a subsequent publication. 



2 Total cross sections 



As already noted, there is a substantial existing literature on vector boson pair production 
in hadronic collisions. As a cross-check of our results, we will compare the values of the total 
cross section obtained using our Monte Carlo (MCFM) with those of Frixione, Nason, Mele 
and Ridolfi 0, Bl ffl. Since the earlier predictions were made at centre-of-mass energies of 
yfs = 16 TeV and ^Js = 40 TeV, we will also provide up-to-date values for pp collisions at 
2 TeV and pp collisions at 14 TeV. These Run II Tevatron and LHC predictions contain the 
latest parton distribution sets |T0, [Tl| as well as more recent electroweak input. 



2.1 Comparison with existing results 

We will first present a comparison with older results on the total cross section. We use the 
structure function set HMRSB |T2"fl , which is common amongst the calculations of [3|, ^]. 
This set corresponds to a four-flavour value of AH? D = 190 MeV. For the purposes of 
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y/s = 16 TeV 
(PP) 


W + W~ 




ZW~ 


ZZ 


MCFM 


Ref. 1 


MCFM 


Ref. [1 


MCFM 


Ref. § 


MCFM 


Ref. § 


Born [pb] 


64.10 


64.11 


14.64 


14.61 


10.26 


10.25 


9.76 


9.75 


Full [pb] 


99.03 


99.03 


27.18 


27.11 


18.95 


18.91 


13.3 


13.2 



Table 1: Total cross section for the various di-boson processes at y/s = 16 TeV. For the 
cases of W ± Z and W + W~ production the electroweak parameters are, m\y = 80.0 GeV, 
m z = 91.17 GeV, = 128 and cos 6^ = m w /mz^, 0- For Z pairs the input is instead, 
m z = 91.18 GeV, sin 2 W = 0.228 and = 1280]. 



comparison we use the standard two loop expression [L3 

a s (Q 2 



Mn(Q 2 /A 2 ) 



6'lnln(Q 2 /A 2 ) 



b ln(Q 2 /A 5 



(1) 



and match to five flavors at m b = 5 GeV, a^m;,) = af^im^). This yields a strong coupling 
at the Z-mass of, 

a s {m z ) = 0.10796. (2) 

The renormalization and factorization scales are chosen to be equal to the average mass of 
the produced boson pair (namely, mw, m z and (mw + m z )/2). 

Finally, since the early calculations present results for the production of two on-shell 
vector bosons, we need to ensure that in our Monte Carlo (in which we produce 4 final-state 
leptons) we use the narrow-width approximation where only doubly-resonant diagrams are 
included. The on-shell boson cross sections can be obtained by dividing out by the relevant 
branching ratios. A full discussion of the other diagrams included in our approach is given 
in Section 0. 

The comparison of the results for a pp collider at centre of mass energies of yfs = 16 TeV 
and yfs = 40 TeV is shown in Tables |l| and § respectively. For the cases of WW and ZZ 
pair production, the values for comparison are taken directly from and M. For ZW ± , the 
results are slightly different from those published in since the contribution from processes 
of the type g + b — > Z + W~ + t with the top quark taken massless (which it is no longer 
appropriate to include) have been removed[]. The comparison between our results and the 
results of refs. 0-[[§] is satisfactory. Apart from their role as a check of our programs the 
results in Tables [I] and || should be considered obsolete. 



We are grateful to P. Nason for providing us with a modified code for these cases. 
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y/s = 40 TeV 
(PP) 


W + W~ 


ZW+ 


ZW~ 


ZZ 


MCFM 


Ref. 1 


MCFM 


Ref. [1 


MCFM 


Ref. § 


MCFM 


Ref. § 


Born [pb] 


148.9 


149.0 


33.56 


33.49 


25.58 


25.54 


23.5 


23.5 


Full [pb] 


254.1 


254.1 


71.86 


71.71 


54.78 


54.67 


34.1 


33.9 



Table 2: Total cross section for the various di-boson processes at ^/s = 40 TeV. Input 
parameters are given in the caption to Table [l|. 



y/s = 2 TeV 
(pP) 


w+w- 


ZW+ or ZW~ 


ZZ 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


Born [pb] 


10.0 


10.3 


1.46 


1.49 


1.22 


1.25 


Full [pb] 


13.0 


13.5 


1.95 


2.01 


1.56 


1.60 



Table 3: Total cross section for the various di-boson processes for the Tevatron Run II. Input 
parameters are given in the text. 



2.2 Results for Tevatron Run II and the LHC 



In order to update the currently available predictions, we use modern values of the particle 
masses and widths, as given in JT3| , 



M w = 80.41 GeV, Y w = 2.06 GeV, 

M z = 91.187 GeV, T z = 2.49 GeV, (3) 
together with a~l^ = 128.89 and present the results using the most recent distributions of 



two popular sets of structure functions, MRS98 [10] and CTEQ5 [fLl] 



The total cross-sections expected at the Tevatron and the LHC are shown in Tables [| 
and |j. We have used a central gluon and as(M z ) = 0.1175 for the MRS98 parton distri- 
bution set (ft08a), whilst CTEQ5M has as{M z ) = 0.118. As before, the factorization and 
renormalization scale fi is set equal to the average of the produced vector boson masses. 
Note that because of changes in the structure functions and a s the modern cross sections at 
y/s = 14 TeV lie well above the old values at y/s = 16 TeV given in Table |l|. 

One can see that in Run II at the Tevatron, the K-factor (the ratio of the full next-to- 
leading order result to the Born level prediction) is approximately 1.3 for each case, whilst 
at the LHC it varies between 1.3 for Z-pairs and 1.7 for ZW~ . The differences between the 
two choices of parton distributions considered in this paper are of the order of 3% in Run 
II, but about 6% at the LHC. 

Fig. [I] shows the scale dependence of the cross section at ^/s = 2 TeV both in leading and 
next-to-leading order using the MRS98 distribution. The growth of the cross sections with 
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= 14 TeV 

(pp) 


W + W~ 


ZW+ 


ZW~ 


ZZ 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


Born [pb] 


81.8 


86.7 


18.6 


19.9 


11.7 


12.5 


12.2 


12.9 


Full [pb] 


120.6 


127.8 


31.9 


34.0 


20.2 


21.4 


16.3 


17.2 



Table 4: Total cross section for the various di-boson processes for the LHC Input parameters 
are given in the text. 

energy is shown in Fig. |2], emphasizing that at high energy vector boson pair production is 
dominated by production off sea partons. Note however that it is still true that a{W + Z) > 
a(W~Z) at the energy of the LHC. 



3 Beyond the zero-width approximation 

Part of the reason for re-evaluating the vector boson pair production cross sections is to 
estimate their importance as backgrounds for new physics processes. In this context the 
tails of the Breit-Wigner distributions may be important. We are therefore motivated to 
go beyond the zero width approximation. We consider all standard model contributions to 
four lepton production, rather than just those proceeding through the production of a pair 
of vector bosons. 

In the zero width approximation q 2 = My, the doubly-resonant diagrams form a gauge 
invariant set. If we wish to move beyond the zero-width approximation, so that we have 
q 2 ^ My, gauge invariance requires that we include all diagrams which contribute to a given 
final state. This problem has been extensively studied in the e + e~ environment where similar 
diagrams contribute to W^-pair production . 



In practice this means that in addition to the diagrams containing two resonant propa- 
gators, we must calculate diagrams containing only a single resonant propagator. For the 
case of IV-pair production, examples of such doubly and singly resonant diagrams are shown 
in Figure [5| To illustrate the gauge-dependence of the individual sets of diagrams one may 
work in an axial gauge, taking care to include both the mixed propagators and additional 
vertices that arise in these gauges. 

Even when we include all the diagrams, a second problem arises when we introduce a 
width for the propagators to avoid on-shell poles. With the modification, 

1 1 



q 2 -M 2 q 2 — M 2 + iMV 

for each of the propagators, the amplitude is no longer gauge invariant because we now have 
a mix of singly and doubly resonant diagrams. The Breit-Wigner form of the propagator 
sums self-energy diagrams which are not separately gauge invariant. Since the resummation 
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Figure 1: Scale dependence of vector boson pair production cross sections. 
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Figure 2: Energy dependence of vector boson pair production cross sections. The scale /i is 
taken to be the average vector boson mass. 



6 




7 



of all diagrams which contribute to a given process is not practical, several models have been 
proposed which allow the introduction of a finite width but preserve gauge invariance. 



For a review of the models, the 'pole-scheme' 

0- 



15 1, the 'fermion-loop scheme' and the 



'pinch technique' 



see for example Ref. 



Here we will adopt the simple prescription 
whereby we use l/{q 2 — M 2 ) for each propagator initially and then multiply the whole 
amplitude by, 



n 

props 



T 



M 2 



q 2 -M 2 + iMY 



which clearly maintains gauge invariance [JT8J . This is the correct treatment for the doubly- 
resonant piece, but mistreats the singly-resonant diagrams, primarily in the region q 2 ~ M 2 
where the doubly resonant diagrams dominate. This is the method that we will use to 
produce our Monte Carlo results in the remainder of this paper. Since the introduction of 
the width represents an all-orders resummation of an partial set of diagrams, there is no 
unique way to include it at a given order. 



3.1 Cross-sections to leptonic final states 

We now present our results for the various leptonic final states using the prescription given 
above. In practice, the observed cross sections are limited by the acceptance of the detectors 
in rapidity and in transverse momentum. We quote cross-sections for particular channels 
with cuts appropriate for Run II and the LHC Specifically, we apply cuts on the transverse 
momentum and rapidity of each lepton, 



p T > 20 GeV , < rj- 



max • 



J 2 Run II 

and a cut on the total missing transverse energy, 

$ T > 25 GeV. (5) 

This suffices for WW production; for the other cases we perform a mass cut 75 < Mi+i- < 
105 GeV. For a more detailed study, one might tailor the cuts individually for each process 
and include detector resolution effects. 

Using the same input as the previous section, but now with the inclusion of the singly- 
resonant diagrams, our full next-to-leading order results with these cuts are summarized in 
Tables H and H. 



4 Examples 

In this section we present two examples to demonstrate the use of the MCFM Monte Carlo. 



8 



Run II 


W~W+ -» e"Pe+^ 


ZW ± -» e-e+Pe* 


ZZ — ► e e + /i, yU + 


ZZ -» e-e+(z^P x 3) 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


« [fb] 


70.9 


73.5 


2.89 


2.99 


1.64 


1.70 


10.7 


11.0 



Table 5: Cross sections (in fb) for various channels at the Tevatron Run II with the cuts 
of (D and (|). 



LHC 


W~W + -> e've+v 


ZW + -> e-e+z/e+ 


ZW~ e"e + i>e- 


ZZ — > e e + /i /i + 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


MRS98 


CTEQ5 


<t? [fb] 


514 


549 


34.0 


36.2 


22.9 


24.6 


12.0 


12.7 



Table 6: Cross sections (in fb) for various channels at the LHC with the cuts of 
4.1 Tri-lepton production at Run II 



and 



One of the 'gold-plated' supersymmetry discovery modes at Run II is gaugino pair-production 
resulting in a tri-lepton signal, 



pp 



(6) 



where Xi is t ne lightest supersymmetric particle. In order to obtain a clean signal, it is 
imperative to have a good understanding of the Standard Model background, which is pre- 
dominantly from leptonic decay of a WZ pair. There have been many previous studies of this 
background (see |3Uj and references therein), primarily performed using the IS A JET [21] and 
PYTHIA event generators. Whilst event generators can have advantages over a fixed 
order calculation, as presently written ISAJET and PYTHIA do not include the 7* contri- 
bution or the interference between the photon and the Z in WZ/j* production. 



In the recent paper P0|| , a set of cuts was proposed to further isolate the trilepton signal. 
However, this study was based on a PYTHIA analysis which includes only the WZ process 
and not Wj* process in assessing the standard model four lepton background. 



Following p0| , we apply the cuts, 



central lepton : 
remaining leptons : 



> 



p T > 11 GeV, \rf\ < 1, 
p T > 7 GeV and p T > 5 GeV, 
25 GeV, 



(7) 



and examine the invariant mass distribution of opposite-sign, same-flavour lepton pairs. This 
differential cross-section (using the MRS98 structure functions at ^/s = 2 TeV) is shown in 
Figure |[ 
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25 50 75 100 125 150 

ixii+i- [GeV] 

Figure 4: The lepton invariant mass with the full Z/^* interference (single hatched) and with 
the Z only (double hatched). Each generated event is binned once for each opposite-sign, 
same-flavour lepton pair. Events which give rise to two entries are binned with half the event 
weight. The total cross-sections in the plot are 49.4 fb and 42.6 fb respectively. 
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We see that as the invariant mass of deviates from the Z-peak, the contribution from the 
off-shell photon dominates. For the case presented in ref. |20| we have m x 122 GeV which 
gives a signal region 10 < M\+\- < 60 GeV. Studies using PYTHIA therefore underestimate 
the standard model background. 



4.2 Higgs production via gluon-gluon fusion 

At the Tevatron the dominant production mode for a (Standard Model or similar) Higgs 
boson is via gluon fusion, gg — > H via heavy quark loops. A natural decay mode for 
140 < M H < 180 GeV is then H — > W^W^ — > leptons/jets, which has been extensively 
discussed in the literature |23], 24]. The most recent of these studies |p4]| , performed using 



PYTHIA, optimized a set of cuts to suppress the SM backgrounds for the di-lepton plus 
missing energy channel and for the like-sign lepton plus jets channel. 

Here we perform a parton-level analysis for the di-lepton + missing energy signal. The 



signal is calculated using the heavy-top effective ggH vertex |25[ and we have applied cuts 



(10)-(16) of which we now describe. In addition to a standard set of cuts, 



p T (e) > 10 GeV, \r) e \ < 1.5, 

p T (/ii) > 10 GeV, p T (/i 2 ) > 5 GeV, \r] \ < 1.5, 

m(ee) > 10 GeV, AR(£j) > 0.4, $ T > 10 GeV, (8) 

there are further cuts to reduce the various background processes. First we apply, 

<j){££) < 160°, 9(££) < 160°. (9) 

where <p(££) is the azimuthal angle in the transverse plane and 8(££) the three-dimensional 
opening-angle between the two leptons. We also impose, 

p T (££) > 20 GeV, cos^ j < 0.5, M T (£^1 T ) > 20 GeV, (10) 

where is the relative angle between the lepton pair transverse momentum and the 

missing transverse momentum and the two-body transverse-mass is defined for each lepton 
and the missing energy as, 

M 2 T {£$ T ) = 2p T {£)$ T {\ - cos^). (11) 
Further di-lepton mass cuts are, 

m(££) < 78 GeV for e + e-,/i + ^, 

m(££) < 110 GeV for e/i, (12) 
and we also cut on the Dittmar-Dreiner angle Q\ for each lepton £ 1; 

-0.3 < cosfl; < 0.8. (13) 
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Finally, we introduce a jet veto, 



veto if > 95 GeV, |^| < 3, 

veto if p% > 50 GeV, |^.| < 3, (14) 
and reject events where either jet is 5-tagged with an efficiency, 

e b = 1.1 x 57% tanh(-^-). (15) 

3o. 1)5 

The results of our analysis (cr cuts ) at y/s = 2 TeV for 140 < m# < 190 are shown in Ta- 
ble |7], where we have employed the structure function set MRS98. Also shown in Table [7| are 
the total cross sections for the various channels which serve as normalizations of our results. 
The primary background in this channel is from W + W~ with smaller backgrounds from top 
pair production?], WZ and ZZ production. For the signal, we choose the renormalization 
scale fj, = rriH] for the di-boson backgrounds we again use the mean boson mass. In the 
case of the tt background we set /i = 100 GeV; order corrections to the total top pair 
production cross-section are small at this scale p6[. 



Our analysis confirms that the WW process is the principal background for Higgs pro- 
duction in the region 140 < inn < 190. Note that we have not included the decay 
W{— > tv — > / + X) in our calculations. The comparison with ref. |24| (where these ef- 
fects are included) is therefore not exact. The effects of r decays are stated to be small in 
ref. [El. 



We find that all the backgrounds due to the di-boson processes are larger than in ref. f24 
primarily because we have normalized to the 0(a s ) cross-sections which are about 30% bigger 
than the Born cross-sections, (c.f. Table |3|). Note however that our signal cross sections are 



also larger than Ref. pl| by about 20%, so the net effect on S/ y B may be small. In addition, 



we find that the background from the P^Z-class of events is about twice as big as ref. [E4 



because of our inclusion of the Wj* contribution, which was left out in ref. ||24|| . Our estimate 
of the tt background is smaller than in Ref. ||24|| , perhaps because we do not include leptons 
from 6-decay. 



5 Conclusions 



We have presented a Monte Carlo program for vector-boson pair production at hadron 
colliders, including for the first time the complete 0(a s ) corrections with leptonic decay 
correlations. 

We have employed this program to calculate the total di-boson cross-sections, first as a 
cross-check with existing results in the literature and secondly in order to provide an update 
of predictions for Run II at the Tevatron and for the LHC, including the latest structure 

2 This process is implemented only at leading order in MCFM. 



12 



Signal (ttih) 


140 


150 


160 


170 


180 


190 


Ocuts (fb) 

0"total (Pb) 


4.36 
0.181 


5.32 
0.206 


6.12 
0.215 


5.15 
0.180 


3.90 
0.143 


2.47 
0.0976 


Background 


WW 


tt 


ZZ 


WZ 


Ocuts (fb) 
^totai (pb) 


185 
13.0 


9.55 
6.82 


2.48 
1.56 


10.4 
3.96 



Table 7: Signal and background cross-sections for a Higgs search with di-lepton final states. 
The t + t~ background is negligible with these cuts. For comparison, the total cross sections 
(without leptonic decays from the vector bosons) are also shown. 

functions and strong coupling. The cross sections are larger than previous estimates in the 
literature. 

The advantage of this Monte Carlo is only realized when cuts are applied to the final-state 
leptons. This is primarily of importance when estimating Standard Model backgrounds to 
new physics, which is especially crucial at the Tevatron. For this reason we have provided 
two examples of such uses in Run II, tri-lepton production as a SUSY signal and a di-lepton 
analysis for an intermediate-mass Higgs search. 



A Appendix: Amplitudes for four fermion processes 

We first introduce a separation of the total n-particle tree-level amplitude into its doubly- 
and singly-resonant components, 

AT = C coup (A^k + <sr) , (16) 

where C coup is an overall coupling factor depending on the di-boson process under consider- 
ation. The contributions -A^dr were first calculated in Jl9| although we will closely follow 
the more recent notation and approach of In particular, the amplitudes are presented in 
terms of particle momenta Pi that are all outgoing, so that momentum conservation implies 
EiPi = o. 

A.l W+W final states 

For this process we label the particles as, 

- M 2 W- (- £ 3 u 4 ) W+ (- £ B v 6 ) , 

where the leptons are not necessarily of the same flavour, and we write the process in this 
manner to remind the reader that q\ represents an outgoing quark (so that when we cross 
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to obtain the desired result, it becomes an incoming anti-quark). Then from |T| we see that 
the non-vanishing doubly-resonant helicity amplitudes for up-quark annihilation are (with 
labels suppressed where possible), 

^drK, «£) = (AT e ' a (h 2, 3, 4, 5, 6) + C^jg^Q., 2, 3, 4, 5, 6)) V U V 56 , 
^Sdr«,«2) = C R , u AT' b (2, 1,3, 4, 5, 6)^56, (17) 

with the sub-amplitudes Ag ee given by equations (2.8) and (2.9) of [IJ]. The T'y are propa- 
gator factors given by, 

Vl ' = T^fe ' (18) 

with M12 = M z and M 34 = M 56 = M w . The couplings that appear in these amplitudes are, 

C conp = (-4^-) , (19) 
\ sin W / 



C i>0 = ±2g rd} sin 2 ^ + (l T 2Q { « } sin 2 ^)P 12 , (20) 
C Rr} = ±2g { . } sin 2 ^(l-P 12 ), (21) 



where Qi is the electric charge in units of the positron charge. For the additional singly- 
resonant diagrams we find, 

AlJ R (u^u*) = 2sm 2 6 w x{ 

(V 56 AT e ' a (3, 4, 6, 5, 2, 1) + V u AT e ' a (Q, 5, 1, 2, 4, 3)) v L>n v L , u P 12 
+ (v 5e AT e ' a (3, 4, 1, 2, 5, 6) + P 3 44 ree,a (6, 5, 3, 4, 2, 1)) (Q u Q e + v L , e v L , u V 12 )} , 
AlJ R (u R , u L ) = X rc sR(^, u r ) (L <-> J2, 1 <- 2) , (22) 

where we have introduced a further set of scaled couplings |]J , 

-1 - 2Q e sin 2 6'vy 2Q e sin 2 6V 

= , v R>e = — , (23) 

sin 2u w sin 2U W 

±l-2Q q sm 2 6 w 2Q q sm 2 e w . 

Vl « = shi2^ ' VR - = ~ sin 2% ' (24) 

vl >» = ^W- (25) 

Note that here the sub-amplitudes needed for the singly-resonant diagrams are exactly those 
introduced in |1] to describe the doubly-resonant diagrams. This is because the diagrams are 
topologically equivalent and (modulo couplings) we can obtain the singly resonant amplitudes 
by simply re-labelling the external legs. The down-quark amplitudes may be obtained simply 
by symmetry, 

<drM = ^ R (ti, ti) (« «-> d, 3 <-> 6, 4 <-> 5) , 

AlJ R (d,d) = AlJ R (u,u)(u~d). (26) 
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As described in |lj , the doubly- resonant amplitudes for the process with an additional gluon 
radiated from the quark line are exactly analogous to fll7D and require the introduction of 
functions .A^tree- However, unlike the 6-particle amplitudes, these functions are not sufficient 
to describe the singly-resonant diagrams. In this case, initial state gluon radiation in the 
singly-resonant diagrams would correspond to final-state radiation in the equivalent doubly- 
resonant diagrams and thus we need to introduce a new set of sub-amplitudes. For a positive 
helicity gluon of momentum pj we find, 

A t ^ K (utulgl) = 2sm 2 e w x{ 

{V m B^ a {l, 2, 3, 4, 5, 6, 7) + V M B? ee > b (2, 1, 3, 4, 5, 6, 7)) v Ltn v Li „V m 

+ {v m B^ h {2, 1, 6, 5, 4, 3, 7) + V u B? cc ' a (l, 2, 6, 5, 4, 3, 7)) (Q u Q e + v L , e v L , u V 127 )}, 

(27) 

where the new functions are defined by, 

B tree,a » (36) (1 1 (2 + 7) |4) (1 1 (3 + 6) 1 5) 

7 ~ (17)(27) S56 t 127 t 3 56 ' 1 ' 

and, 

B tr Ce , b -i (26) [45] [(12) (3|(5 + 4)|1) - (27) (3|(5 + 4)|7)] 

(17)(27)s 3 4ti2 7 t 345 

Our definition of the spinor products follows ref. |TJ]. The remaining amplitudes are obtained 
as above, with the additional (negative helicity) gluon corresponding, as in the doubly- 
resonant case, to the operation — flip 1 defined in [jy]. As will be the case for all the amplitudes 
in this appendix, the loop contributions for the singly-resonant diagrams may be simply 
obtained by following the prescription (3.20) of [|TJ. 



A.2 W ± Z final states 

Here we label the two processes slightly differently, 

-> giferhto)^^^), 

-> qitaW + (-+u 3 h) Z (30) 
in order to simplify the form of the amplitudes and the overall coupling (c.f. ref. |16|) is, 

/ „2 \ 2 



c 



coup 



\sin 9 W ) 

The doubly-resonant amplitude for a left-handed decay of the Z is, 

^DR(gi,?2,4) = 4 ree ' a (i,2,3,4,5,6)K <?2 ^, e p 56 + g g2 g e )p 3 4 

+ 4 ree '°(l, 2, 6, 5, 4, 3) (v L>qi v L>e V 5e + Q qi Q e ) V u 

± 4 ree ' 6 (l, 2, 3, 4, 5, 6) {v L>e cot 6 W V 5G + Q e ) V 12 V 3i , (31) 
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where the masses in the propagators are now M\i = M34 = M\y and M§§ = Mz- For W + 
production we have q\ = d, q 2 — u and the third line has a positive contribution, whilst W~ 
corresponds to q\ = u, q 2 = d and a negative sign. This reduces to the form given in [|]] if we 
set Q e = to neglect the virtual photon diagrams. The right-handed amplitude is obtained 
by a symmetry transformation, 



A T Sr(QU Q2, O (VL,e <"> V R>e , 5 ~ 6) . 



(32) 



The singly resonant diagrams are somewhat more complicated. With a Z propagator we 
can couple both electrons and neutrinos, while 7* may only couple directly to the electrons. 
In addition, if the final state electrons are both left-handed, there is a contribution from a 
diagram containing two W propagators. In total we obtain, 



AZZ(qx,q2,£%) 



= (^ L)1 4 ree ' a (3, 4, 6, 5, 2, 1) + v L)2 A%°*' a (3, 4, 1, 2, 5, 6)) V 5& Vi 2 v L>e 
+ ( Cl 4 ree ' a (3, 4, 6, 5, 2, 1) + c 2 4 ree ' a (3, 4, 1, 2, 5, 6)) V 12 Q 2 e 

+ (c^r ' a (6, 5, 1, 2, 4, 3) + c 2 Af e ' a (6, 5, 3, 4, 2, 1)) (33) 



where the newly introduced couplings vl,i and Cj depend upon the process under considera- 
tion and are given by, 

W + : ci = 0, c 2 = 1, = w L>n , w Li2 = v L>e , 

W~ : Ci = 1, C 2 = 0, f L ,l = V Lj e, f L , 2 = fi,n- (34) 

The right-handed contribution is similar, but does not contain the corresponding final term, 
Ajnidi, 92, = (v LA AT e ' a (3, 4, 5, 6, 2, 1) + ^, 2 < ee ' a (3, 4, 1, 2, 6, 5)) V 56 V 12 v R , e 

+ ( Cl < ee * a (3, 4, 5, 6, 2, 1) + c 2 4 ree ' a (3, 4, 1, 2, 6, 5)) V 12 Ql (35) 

We now turn to the graphs including gluon radiation. For the doubly resonant contribution 
with a positive helicity gluon we find a similar structure, 

A^ n (qi, q 2 , 4 gi) = Af ee ' a (l, 2, 3, 4, 5, 6, 7) {v L>q2 v L>e V 56 + Q q2 Q e ) V u 

+ 4 ree '°(l, 2, 6, 5, 4, 3, 7) (v L , qi v L>e V 56 + Q qi Q e ) V 34 



± A 



tree,b, 
7 1 



1, 2, 3, 4, 5, 6, 7) (v L>e cot fl^e + Qe) V l27 V u , (36) 



whilst the singly resonant pieces again require the new amplitudes, 



A^ R (q 1} q 2 ,£^gi 



+ 



+ 



v B^il, 2, 3, 4, 5, 6, 7) + v L , 2 B? ee ' b (2, 1, 6, 5, 4, 3, 7) ) V 56 V 127 v L 



Cl B^ a (l, 2, 3, 4, 5, 6, 7) + c 2 B* ee ' b (2, 1, 6, 5, 4, 3, 7) ) V 127 Q 



ci5f e ' fe (2, 1, 3, 4, 5, 6, 7) + c 2 ^ ree * a (l, 2, 6, 5, 4, 3, 7) 



tree. a / 



^127^34_ 

2 sin 2 0- 



,(37) 



w 



VL lB ^c,a^ ^ 3j 4) 6? 5) ?) + ^ >2jB tree, 6(2) ^ ^ ^ ^ g^ y) j 
CijB tree,a (1) ^ g^ ^ ^ ^ ?) + ^tree.b^ ^ ^ ^ ^ g^ ?) \ (gg 
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The remaining amplitudes, with the helicity of the gluon reversed, can be obtained by the 
transformation, 

& 4, <?f) = -Aj R ( gi , q 2 , £ 6 , g+) (4 <- 6, 3 <- 5, (oft) - H,^™'" - • 

(39) 
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